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1. Introduction 

The early phase of the Epoch of Reionization (EoR) or Cosmic Dawn (CD) is, by definition, 
believed to have started with the formation of the first stars, even though the EoR may also have 
been fuelled initially by somewhat unconventional radiation sources such as annihilating dark mat¬ 
ter clumps (e.g. Spolyar et al. 2009). The first stars are believed to have formed within minihalos 
(halos with T wir < 10 4 K, or 10 4 < M/Mq < I0 7 8 ) and to have been massive, emitting a sufficient 
amount of UV radiation to ionize the surrounding hydrogen and helium. These stars were born at 
zero metallicity, with their formation predominantly regulated by I E cooling, and thus they are also 
Population III (Pop III) stars. Until recently the conventional view, from high-resolution numerical 
simulations, was of massive (> 1OOM 0 ) stars forming in isolation (e.g. Abel et al. 2002; Bromm 
et al. 2002; Yoshida et al. 2006): one Pop III star per minihalo. This paradigm has undergone 
recent revision due to newer simulations in which smaller-mass (~ 10 — 40 M 0 ) stars form binary 
systems (Turk et al. 2009; Stacy et al. 2010; Greif et al. 2011a), although the universality of this 
result is unclear (see Hirano et al. 2014 for a wide spectrum of Pop III stars in ~ 100 minihalos). 
At any rate, a change in the initial mass function (IMF) of Pop III stars affects the hardness of 
the spectral energy distribution (SED) of these stars, and thus the CD and EoR modelling as well. 
Stellar binaries may evolve into X-ray binaries that may, if efficient, result in a smoother ionization 
structure of the intergalactic medium (IGM) than that by UV sources, due to a much longer mean 
free path than that of typical UV photons (e.g. Haiman 2011; Mesinger et al. 2013 and references 
therein). Furthermore, the interstellar medium heated by supernova explosions can cool through 
a combination of Bremsstrahlung and metal line cooling, producing soft (<keV) X-rays that can 
efficiently heat the high-redshift IGM (Pacucci et al. 2014). 

The theory of the formation and evolution of the first stars, therefore, is crucial in modelling 
CD and even the EoR, because the metallicity built up after the death of Pop III stars will give way 
to the formation of Population II (Pop II) stars, believed to be the main drivers of reionization. A 
major feedback that affects their formation is the dissociation of EE molecules - cooling agents - 
inside minihalos by Lyman-Werner band radiation (e.g. Haiman et al. 2000; for recent simulations 
with this effect in > 100 Mpc box see e.g. Ahn et al. 2012 and Fialkov et al. 2013). Stars would 
not have formed either in minihalos that were embedded in photoionized gas (T > 10 4 K), since 
the increase in pressure would have suppressed minihalo gas accretion (Jeans-mass filtering). 

There also appeared a noteworthy discovery having to do with the nature of large scale fluctu¬ 
ations. Tseliakhovich & Hirata (2010) used a peak-background split scheme to study the nonlinear 
evolution of the baryon/dark-matter velocity offset, which was seeded at the recombination epoch, 
and found that hydrodynamics at scales relevant to minihalo formation should have been affected 
(see more details of its astrophysical impact in the subchapter by Maio). The typical offset is found 
to be e.g. ~ 1 km/s at z ~ 20, corresponding to Mhaio ~ 1.5 x 1O 5 M 0 , when minihalos were al¬ 
ready abundant. While determination of the halo mass and environmental dependence of this effect 
awaits more detailed study, simulations already show an increase in the threshold mass for star 
formation and a delay in halo collapse (Stacy et al. 2011; Greif et al. 2011b; O'Leary & McQuinn 
2012). More importantly, this nonlinear effect may shock-heat the gas globally, introducing a new 
feature into the power spectrum that corresponds to the relative velocity fluctuations (McQuinn & 
O'Leary 2012). When the global shock-heating is efficient, this new feature could even dominate 
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over the power spectrum of linear density perturbations and amplify the baryon acoustic oscillation 
(BAO) feature. Due to the uncertainties in the efficiency of this shock-heating, however, the actual 
amplitude of the power spectrum due to relative velocity fluctuations is still uncertain. 

From an observational perspective on CD and the EoR, especially in terms of 21-cm power 
spectrum analysis, it is convenient to mark three prominent epochs: (1) the Lya-pumping epoch, 
when the IGM is strongly coupled to 7k through Wouthysen-Field effect with a high Lya intensity, 
(2) the X-ray heating epoch, when the IGM is gradually heated to beyond 7 cmb by X-ray heating 
and (3) the EoR, when H II bubbles in cosmological scales form in a patchy way. It is generally 
believed that CD commences with the Lya-pumping epoch, followed by the X-ray heating epoch, 
and finally occurs the EoR, whose sequence is rather robust unless 1-2 order-of-magnitude changes 
in the fiducial astrophysical parameters are allowed (Furlanetto 2006; McQuinn & O’Leary 2012; 
Mesinger et al. 2014). Very conveniently, the early phase of each epoch boosts the spatial fluctua¬ 
tion of 8Tb, sequentially dominated by the patchiness in the Lya intensity, the IGM temperature, 
and the ionized fraction, respectively, extending the observational window to very high redshifts. 

While challenging, it should be possible to observe individual objects via tomography, which 
may be as small as the first galaxies. When UV sources are embedded in an IGM colder than the 
CMB, a “Lya sphere” forms around the radiation source and a strong 21-cm absorption trough 
forms. When UV sources are accompanied by X-ray sources, the central, heated region is seen in 
emission, which is still surrounded by an absoiption trough or an absoiption plateau (e.g. Tozzi 
et al. 2000; Cen 2006; Chuzhoy et al. 2006; Chen & Miralda-Escude 2008; Alvarez et al. 2010). A 
simulation of the formation and evolution of Pop III and II stars in a rare density peak, together with 
the X-ray binaries, shows that this ubiquitous feature will be observable by SKA2-LOW with 1000- 
hour integration, and marginally by SKA1-LOW (Ahn et al. 2014) with very aggressive integration 
at z < 15 (Section |3|). 

2. Power Spectrum Analysis 

Because many independent modes with similar wavenumbers ( k = |k|) are averaged over to 
form the 3D power spectrum, P(k), power spectrum analysis has higher sensitivity than imaging 
in general. In addition, nature cooperates in such a way that three important physical processes - 
Lya pumping, X-ray heating and patchy reionization - may sequentially boost 21-cm fluctuations 
to give a roughly constant S/N ratio over redshift, even though the foreground increases rapidly 
towards high redshift (e.g. Mesinger et al. 2014). Overall, the signal and the foreground noise grow 
roughly at a similar rate to yield S/N~ 10 ' 5 , although the S/N ratio should eventually become very 
small (~ 10 7 ) at z ~ 30 (Pritchard & Loeb, 2008). 

The very first light from astrophysical sources to which the IGM is exposed is the stellar con¬ 
tinuum below the hydrogen Lyman limit, or photons with energy hv < 13.6eV, because ionizing 
photons (photons above the Lyman limit) are quickly absorbed by the neutral hydrogen and helium 
around while continuum photons are not. However, there is always some fraction of this contin¬ 
uum that is absorbed in Lyman series resonances as photons redshift over cosmological distances 
(e.g. Ahn et al. 2009). After absorption, cascades over the energy levels and multiple scattering of 
reemitted photons will convert about ~ 30% of those absorbed photons to Lyman-a photons (Hi- 
rata 2006; Pritchard & Furlanetto 2006). Even in the presence of density inhomogeneities, a point 
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Figure 1: (top subpanels) Evolution of D 2 (k) at k = 0.1 /Mpc and sensitivity of various planned observations 
including SKA1-LOW (We note that the apparent SKA1-LOW sensitivity cut at 20 is due to the particular 
choice of instantaneous bandwidth, A z = 0.5, which limits the number of available k = 0.1/Mpc modes; we 
stress that bandwidth choices are arbitrary and that SKAl-LOW’s planned bandpass extends out to z = 28.). 
Three peaks appear when fluctuations in 57/ are boosted due to the efficient Lya-pumping, the X-ray heating 
and the patchy reionization, subsequently from high to low redshifts (Mesinger et al. 2014). 


source produces a nearly spherical profile (Vonlanthen et al. 2011) that shows a step-wise feature 
decreasing more steeply than 1/r 2 (Pritchard & Furlanetto 2006). As usually assumed, the IGM 
is still much colder than the CMB, and thus the regions where the Lya pumping becomes strong 
(N a > N a ,th — 10 10 (yyr) s~ 'em 2 Hz 'sr'" 1 ) will show strong absorption, or 57/ ~ — lOOmK. 
Source clustering and the steep radial profile of the Lya intensity around each source are the dom¬ 
inant sources of fluctuations in 57/ at earliest times. Soon thereafter, saturation occurs, and fluc¬ 
tuations decrease as sources grow in abundance and N a 2> 1 everywhere (Fig. Q]). Changes in the 
spectral energy distribution (SED) of sources can affect the resulting power spectrum as demon¬ 
strated by Santos et al. (2011) (Fig. {2|), while other systematic effects such as the star formation 
efficiency might make the prediction somewhat degenerate. 

In many models, the Lya pumping epoch is followed by the X-ray heating epoch. X-ray 
photons have a much longer mean free path than UV photons, and partially ionize the IGM by 
leaving behind energetic electrons as they traverse cosmological distances. These electrons can 
then ionize, excite, and heat the IGM. During the early phase of the X-ray heating epoch, the overall 
ionization level caused by X-rays is usually very small unless the X-ray emissivity is very high (see 
e.g. Mesinger et al. 2013). Thus the dominant contribution from X-rays to fluctuations in 8Tb is 
inhomogeneous heating, which will boost the power spectrum again (Fig. [I]). Another contribution 
to fluctuations is due to inhomogeneous Lya pumping from excited hydrogen atoms (see e.g. Chen 
& Miralda-Escude 2008), though its contribution is usually subdominant compared to the Lya 
pumping induced by stellar UV photons (but see e.g. Ahn et al. 2014 for a significant contribution 
when X-ray source formation is very efficient). It is usually expected that X-ray heating epoch 
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Figure 2: 21-cm temperature power spectrum at z — 20.3 for a few Ly a emission models, based on the 
simulation results by Santos et al. (2011). Power-law SEDs are used: the number of photons per frequency 
emitted per stellar baryon e v =Av a , where A is tuned to produce the assumed integrated number of photons 
between Lya (10.2 eV) and the Lyman limit (13.6 eV), or Ny (legends are self-explanatory), a = —0.9 
and a = 0.29 roughly represent Pop-II and Pop-Ill type SEDs, respectively. Cases with identical Ny (black 
solid and black dashed curves) are almost degenerate. Pushing the minimum halo mass hosting stars down to 
M m j n = 10 6 Mq, with a = 0.25 and Ny = 5000, renders the overall shape of the power spectrum different from 
the other three cases, which all have M m ; n = 10 s /VC . Error bars in blue corresponds to the full SKA2-LOW 
collecting area, and the ones in red to the 10% of the collecting area (Santos et al. 2011). 


produces the highest peak in the large-scale power spectrum in its redshift evolution (Fig. [TJ). Later, 
X-ray heating will become very efficient everywhere in the Universe such that 8Tb > 0 and 7)gm 
7cmb but with the ionization state still low, once again saturating the fluctuations and resulting in a 
significant decrease in the amplitude of the power spectrum (Fig. |I|). This rise, peak and fall in the 
power spectrum amplitude results in the second “bump” seen in figure [1} Note that this provides the 
best window for high-z cosmology, because 8 Tb becomes proportional to the baryon density and 
nearly insensitive to other fluctuations. In addition, separating the matter power spectrum from that 
of the radiation fields, utilizing the impact by the peculiar motion (or the so-called jU-decomposition 
scheme which is very similar to that used in galaxy survey programs), may be possible (McQuinn 
et al. 2006; Mao et al. 2012; see also the subchapter by Pritchard). 

Finally comes the main EoR, in most models. The power spectrum is boosted again, due to 
the patchiness of H II regions where 57), = 0 and neutral regions where 8 Tb is proportional to the 
baryon density. Final saturation and power spectrum suppression occurs as bubbles overlap and 
the EoR approaches its final phase (e.g. Iliev et al. 2012), to generate the final bump in the power 
spectrum evolution (Fig. |l]). Of course, one should allow more model variants. Currently, the X- 
ray emissivity at CD and during the EoR is still largely unknown, other than the weak constraint 
coming from the soft X-ray background (e.g. Dijkstra et al. 2004; McQuinn 2012), the z ~ 2.5 QSO 
absorption lines and the reionization histories of HI and He II (McQuinn 2012). This allows a wide 
model variation in the X-ray heating epoch, and thus the 3-bump feature in the power spectrum 
(Fig. [T]) is not ubiquitous among reionization models if extreme cases are allowed (Pritchard & 
Loeb 2012). The halo mass spectrum responsible for reionization is not well constrained either, 
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Figure 3: (a:left) Distinct ring structure due to the step-wise Ly a pumping around a point source. - 5'/}, x r 2 
at z = 13.42 is plotted with an arbitrary unit, to show these distinct rings (pointed by arrows) as boundaries 
of sudden change in color: each annulus bound by adjacent rings, where Lya flux decreases as 1/r 2 , shows 
about the same color. The box size is 137 Mpc coming and the angular size is 51' (Vonlanthen et al. 2011). 
(b:right) Radial profile of an isolated radiation source of the UV and the X-ray. From inside out, 81), = 0 (H 
II region), 8Tb > 0 (X-ray heated region), 87), < 0 (Lya-pumped region). 

allowing models with a significant contribution from minihalos (e.g. Wyithe & Cen 2007; Ahn 
et al. 2012). 

The sweet spot for power spectrum analysis exists around k ~ 0.1 /Mpc, guaranteeing a large 
S/N over a wide range of redshifts while being mostly unaffected by the filtering of low-/: fore¬ 
ground modes. SKA 1-LOW (and the 50% capability of SKA 1-LOW as well) has a large enough 
S/N ratio for observing P(k ~ 0.1/Mpc) for z < 20 with 1000-hour integration, and SKA2-LOW 
will allow one to achieve similar S/N ratio for even larger k’ s with more extended redshift cover¬ 
age, z < 28 (Fig. [lj see also Baek et al. 2010 for the evolution of other wavemodes). As is usually 
the case for CD and the EoR, pushing the observable redshift to z ~ 28 or v ~ 50MHz is strongly 
favored, in order to probe their early phases. 

3. Tomography 

As mentioned in Section imaging is more difficult than power spectrum analysis. While 
the power spectrum contains a lot of information on the underlying astrophysics and cosmology, 
one should note that the 3D tomography (imaging) is indeed a more direct observation because 
one cannot invert P(k) to obtain the actual image unless the 21cm field is Gaussian (see also the 
subchapter by Mellema). 

Imaging can probe physical properties of radiation sources. One interesting feature is that 
due to the step-wise decrease of Lya intensity (or Lya pumping) in addition to the geometrical 
dilution, 1/r, the radial profile of the 21-cm signal around an isolated source will also exhibit a 
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Figure 4: Realistic imaging forecast of the early objects with the first stars (Population III stars). Population 
II stars and Population III X-ray binaries, (Ahn et al. 2014). (a:top left) 21cm map (mK) from a simulated 
model of Pop III stars+Pop II stars without X-ray sources in a rare density peak inside a 40 Mpc (comoving) 
box at z= 15, smoothed with © = 2' and Av = 1MHz. The angular size of the box shown is 13' at this 
redshift. (b:top right) 21cm map (mK) from a simulated model of Pop III stars+X-ray binaries+Pop II stars 
in a rare density peak inside the same box as (a) at z = 15, smoothed with 0 = 2' and Av = 1 MHz. (cibottom) 
Imaging sensitivity of the 2-km core SKA2-LOW (Mellema et al. 2013) with 0 = 2' and Av = 1 MHz, for 
integration times of 1000 hours (black, solid curve) and 10000 hours (red, dotted curve), against varying X- 
ray binary model predictions (points). It is difficult to image individual objects by SKA1-LOW unless nature 
cooperates: a few models are marginally observable at z < 15, as seen by the sensitivity of SKA1-LOW with 
the same smoothing filters and a 1000-hour integration time (blue, dashed curve). 


step-wise variation, or multiple “rings” (Vonlanthen et al. 2011; Fig. |3]a). Detecting the multiple¬ 
ring feature from one radiation source remains, however, very difficult even with SKA2-LOW, and 
while stacking many such profiles might allow the detection of the rings, one should know the 
source redshift accurately, because stacking sources from other redshift planes will quickly smear 
out the ring feature. 

A more notable and traceable feature from an isolated source is (unless the IGM is saturated 
7s 3> 7 cmb due to strong X-ray heating) 8T/, varying even more steeply than 1/r 2 . This is due 
to the step-wise variation of Lya pumping rate, resulting in strong absorption nearby and weak 
absorption further away. This “Lya blob” should be ubiquitous around isolated sources before the 
full saturation of Lya pumping occurs everywhere, which is analogous to individual H II regions 
in the patchy reionization process before the full overlap occurs (Fig. [5|t). 
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It becomes even more interesting when the UV source is accompanied by an X-ray source. 
Typical emissivity of X-ray sources associated with UV (stellar) sources is relatively small (com¬ 
pared to proposed Population III X-ray binaries), and the corresponding X-ray heating zone (where 
7s > 7 cmb) is smaller than the Lya blob when the IGM were colder than the CMB before expo¬ 
sure. In addition, the central region will be ionized mostly by UV sources, such that 8Tb is zero 
at the center, positive outside, and negative even further outside (Fig. |b). Note that probing the 
signature of individual objects is likely to be possible in only a relatively narrow range of frequency 
(redshift), since rapid growth in the abundance of galaxies will soon wash out such a signature. 

This can occur even when there is a strongly clustered set of sources: Xu et al. (2014) sim¬ 
ulated the formation of Population III and II stars and X-ray binaries inside a high-density peak 
(“Rarepeak”) in the early Universe, and Ahn et al. (2014) postprocessed the simulated data and 
obtained a 21-cm map which shows a much more extended 8 7), profile than that of a single object 
(Fig. [4a, b). A similar prediction for the spatially extended profile also exists for high-z QSO sys¬ 
tems, very sparsely spaced in the Universe (Tozzi et al. 2000). In case of the Rarepeak, SKA2-LOW 
performing 1000-hour integration of a tracked field with the 0 = 2' and Av = 1MHz smoothing fil¬ 
ter will enable its (their) detection regardless of the model variation, ranging from the all-absorption 
trough case (without X-ray binaries) to the fully X-ray heated case (with efficient X-ray heating) 
with a high S/N ratio (Fig. ^|c). 

As it is easy to estimate the number density of such peaks a priori, one can forecast the odds 
to find those objects inside a dedicated field of observation. For Rarepeak, which is a 3.5 a density 
peak when the linear density field is filtered at R ~ 3Mpc comoving scale, its number density 
is roughly l() 6 Mpc \ A tracking volume of FOV(SKAl-LOW)x(110 — 75)MHz ~ 5° x 5° x 
35MHz (where FOV denotes the field of view), which covers the redshift range z = [ 12 — 17], 
will host about 600 of those high-density peaks. Thus detecting the Lya-blob feature of these 
objects at z < 15 through the SKA2-LOW 1000-hour integration is bound to succeed, while with 
the SKA1-LOW array configuration, about a 16000-hour integration is required to achieve the 
same sensitivity. Alternatively, we may target the few cases which can be marginally detected by 
SKA1-LOW with a 1000-hour integration with the same smoothing scheme (Fig. |4|). While lower- 
redshift detections become much easier, at that time the distinct isolated feature may be erased 
by signals from other, more abundant density peaks if similar emissivity is assumed. In either 
case, since the exact epoch when this occurs is model-dependent, it is necessary to explore a wider 
range of the astrophysical parameter space than has been investigated by Ahn et al. (2014) and still 
consider the possibility to observe individual objects at lower redshifts by SKA1-LOW. It is also 
advisable to carry out imaging with SKA 1-LOW even when it has reached 50% capability of its 
final phase, because again there may exist objects whose size and signal strength are large enough to 
be successfully imaged. For example, AGN+galaxy systems (Tozzi et al., 2000) or even rarer peaks 
than have been studied by Ahn et al. (2014) may appear in such an imaging observation, although 
a relatively low-redshift range, z — 13, may be adequate due to the low sensitivity expected. 

Aside from the high-redshift nature of these early stars/galaxies, the unique feature of the com¬ 
posite H II region + X-ray heated region + Lya blob is not something one expects from the main 
EoR imaging. The main EoR imaging will focus mostly on large H II bubbles (R > 20Mpc) with 
8 Tb = 0 before filtering, and the neutral region will simply fluctuate according to the underlying 
density fluctuation. If so, tomography will tell us about source properties during the peak of EoR 
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(see the subchapter by Mellema and the subchapter by Iliev), and cosmological information during 
the early phase of EoR (see the subchapter by Pritchard). 

4. Conclusion 

We have briefly reviewed theories of high-redshift astrophysics and their observational as¬ 
pects during CD, when the first galaxies formed. The first galaxies significantly affect the IGM, 
leading to signatures that will be detectable through the observation of the 21-cm background. 
With new developments in the theory of first-star formation, the number of feasible models of this 
early epoch is ever increasing. These new developments include the initial mass function of Pop 
III stars, the spectral energy distribution of stars and X-ray sources, the byproduct of stars, and 
dynamical/radiative feedback effects, to name a few. SKA 1-LOW and SKA2-LOW have the capa¬ 
bility of probing first galaxies and their impact on the IGM through power spectrum analysis and 
imaging. 

SKI-LOW will be able to observe the signature of these objects mostly through the 21-cm 
power spectrum, with better sensitivity and thus deeper target-redshift than most SKA precursors. 
Three prominent epochs, which arc the Lya pumping epoch. X-ray heating epoch, and the main 
EoR will be observed as three distinct evolutionary phases in the amplitude of the 21-cm power 
spectrum in the wavenumber range around k ~ 0.1/Mpc. The shape of the power spectrum in a 
wide dynamic range of k might be able to tell us about the properties of the early objects, such as 
their SED. 

SKA2-LOW will be able to carry out imaging, as well as the power spectrum analysis at 
higher sensitivity than SKA1-LOW. Imaging will help to constrain source properties much better 
than power spectrum analysis, because of model degeneracy and loss of information arc inherent 
in the power spectra. Clustered sources of UV and X-ray sources may show their unique signature 
of the composite absorption and emission of 21-cm lines against the CMB, as long as their angular 
scale becomes large enough to achieve high S/N with SKA2-LOW. For serendipity, however, it is 
advised to carry out imaging also with SKA1-LOW. 

There exists tension between high-redshift astrophysics and cosmology. Even at CD, the 21- 
cm signal of an astrophysical origin may swamp that of a cosmological origin. Nevertheless, there 
exists the possibility that the cosmological signal may be boosted significantly, if mechanical en¬ 
ergy conversion occurs on cosmological scales (McQuinn & O’Leary, 2012). In either case, the 
best window for cosmology lies at the highest redshift range planned for SKA (z ~ 28). 
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